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Abstract 

We argue that the superpartner of the Dirac right-handed neutrino is a 
prime candidate for dark matter created from a 'mattergenesis' mechanism. 
We show that due to the smallness of the Yukawa couplings, a right-handed 
sneutrino density created in the early Universe would not be erased by an- 
nihilations, which remain out of thermal equilibrium. It would also not be 
drowned by a later, additional production of right-handed sneutrinos, as the 
relic density of the non-thermal right-handed sneutrinos is found to be gener- 
ally negligible compared to the observed dark matter density. Mild constraints 
on sneutrino masses and trilinear SUSY-breaking couplings are obtained. Pos- 
sible mattergenesis scenarios are also mentioned 



1 Introduction 

While the energy density of baryonic matter, dark matter (DM) and dark 
energy are being measured more precisely, the very nature and origin of dark 
matter remain open questions. In the popular neutralino dark matter sce- 
nario, the dark matter density is equated to the LSP (lightest supersymmetric 
partner of the MSSM) neutralino relic density [1,2]. Indeed, it is often consid- 
ered to be one of the many appealing aspects of the MSSM that it naturally 
provides a dark matter candidate. This scenario, however, fails to address 
the intriguing question of the ratio of the dark and baryonic matter energy 
densities. We now know these to be given by [3] 

n b o.o4 1 1 

The baryon relic density is known to be much smaller than 4% of the crit- 
ical energy density (see eg. [4]) and for this reason it is necessary to find a 
baryogenesis mechanism that can create the right amount of baryons. In the 
common neutralino dark matter scenario, obviously, no such thing is necessary 
for DM since, as we said, the relic density of neutralino is directly the DM 
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density. Since the origins of both types of matter are unrelated, the similarity 
in their present-day densities appears simply as a coincidence. This problem, 
sometimes called the ^DM/^b problem, was first mentioned by [5], and has 
received increasing attention [6-16]. The usual route to tackle this problem is 
to try and create a mechanism that generates both kinds of matter, baryonic 
and dark, at the same time. 

It has been suggested before (eg. in [9,15]) that candidates for mattergenesis- 
induced DM should generally have weak or even super-weak interactions with 
the 'visible' sector. If the candidate never thermalises, the asymmetry cre- 
ated by mattergenesis will not be erased or reprocessed at later times. In this 
case, the smallness of the couplings acts as a built-in protection of the DM 
asymmetry. A condition that other possible sources of DM stay small is also 
added. Another protection mechanism has been suggested in [8], where this 
time the candidate is thermal in the early Universe but freezes-out at some 
temperature T ~ m£)j\,//20, creating a low relic density. The observed DM 
density (and baryon density) is created after freeze-out by the decay of a heav- 
ier particle which couples to both dark and baryonic matter. For this reason 
the mechanism has been called the 'late decay' scenario. In both cases, the 
DM asymmetry is created at a time when the DM candidate is out of thermal 
equilibrium with the plasma, and will remain so. Common to both cases is 
also the fact that other sources of the DM candidate have to be kept small. 

While there is no candidate in the Standard Model or even in the MSSM 
that has a low relic density and never thermalises, once these models are aug- 
mented by right-handed (s)neutrinos to explain the observed neutrino masses 
the situation is somewhat different. We will argue here that the MSSM+Dirac 
right-handed (s)neutrino offers a prime DM candidate for mattergenesis sce- 
narios in the sterile sneutrino LSP. Due to the smallness of its Yukawa cou- 
pling, the right-handed (RH) sneutrino is completely decoupled from the ther- 
mal plasma and has typically a low relic density. In the 'late decay' scenario, 
the condition that the heavy particle should decay after the candidate freezes- 
out would be lifted. This implies that for any mattergenesis scenario creating 
the Dirac RH sneutrino as DM, the DM asymmetry would neither be erased 
by fast annihilation processes nor be drowned by later freeze-out creation. 

In the following we will first introduce a Dirac mass term for the (s)neutrino 
in the MSSM and extract from the model the RH sneutrino interactions. In 
section [3] we will show that under fairly mild constraints on left- and right- 
handed sneutrino mixing the RH sneutrino remains out of equilibrium. We 
will obtain the RH sneutrino relic density numerically in section 2] and show 
that it is much lower than the observed DM density for natural choices of 
parameters. Throughout we shall suppose that a hypothetical mattergenesis 
scenario has created RH sneutrino DM and baryons in quantities mentioned in 
eq.(pQ). Such mattergenesis scenarios already exist in the literature: an Affleck- 
Dine mattergenesis scenario with RH sneutrino DM was proposed in [16] and 
in [17], and earlier an 'early-decay' mattergenesis scenario was suggested in 
[9], where it was noticed that the simplest implementation might be in the 
MSSM+i>r. We will show in the last section that indeed the parameters these 
scenarios necessitate do imply a non-thermal right-handed sneutrino with low 
relic density. 
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2 RH sneutrino interactions 



We add to the MSSM model a RH neutrino superfield N which is given a 
Dirac mass term in the superpotential W: 

W D AI/e^-HV'N , (2) 

where L is the left-handed (LH) lepton doublet superfield and H u is the up- 
type Higgs superfield. We also now have new SUSY-breaking terms 

Vsb = ml L v L *v L + ml R i/ R c *i/ R c + (aXh u W R + /i.e.) (3) 

where a is a mass dimension trilinear coupling, and h u (I) is the up-type higgs 
(lepton) SU(2) doublet. The RH sneutrino has been introduced as a gauge 
singlet and as such will only have a handful of interactions. F-terms are 
the source of 4-point interactions and higgsino exchange (keeping only terms 
involving the RH sneutrino): 



dW 2 x - d 2 W , , 

i,j=L,H u ,N R 3 



i=L,H u ,N R 

= -A 2 (ee*i> R v R * + KK*v R v R * + v~ L v~ L *v R v% + 

+\{D c R H u + e-D R H u \ L ) , (4) 

where all <3?'s (resp. V' s ) stand for the scalar (resp. fermionic) part of the 
superfields, and ht'° is the higgs boson, while Hu'° is the higgsino. The 
trilinear term in the SUSY-breaking Lagrangian provides higgs interactions: 

£sb,v r . = a\ {h+ev c R - h® u v L i> R + h.c.) . (5) 

These are all the interactions the sterile sneutrino has. However once the 
electroweak symmetry is broken, the higgs acquires a vev, and the left- and 
right-handed sneutrino mix0. Indeed, considering this time only mass terms 
in the Lagrangian: 

£mass = -ml R i7 R *v R - ml~ L i/ L *i7 L - a Xv sin $v L v c R - a* Xv sin (3v* L v R (6) 
where the up-type higgs has been replaced by its vev, 

((K),(h° u )) =(0,v sin/3) , (7) 
the usual procedure leads to the following mass-eigenstate RH sneutrino: 



v R u 



aXvsinfi) + fj7i^ — 



[{ m l~ L ~ m l R ) vr - aXv sin^x] . (8) 



: It could also be said that before the electroweak phase transition, left- and right-handed 
sneutrinos mix via 4-point interactions, higgsino exchange or higgs exchange. To call these a 
'mixing' phenomenon or simply interactions is in fact a matter of timescale; would they be fast 
enough, one could say that the RH sneutrino effectively has a left-handed part. We will see however 
in the next section that left-right mixing is in fact very slow compared to the Universe expansion 
before the electroweak phase transition, and indeed never equilibrate. This need not be the case 
once the higgs has acquired a vev; hence our choice of wording. 
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For simplicity reasons, we will take sin (3=1 from now on. As can be expected, 
we will show in the next section that for the RH sneutrino never to thermalise, 
it is necessary that the LH (active) part of it be small compared to the sterile 
part ; in this case eq.([5|) can be rewritten 

\av 

vrm ^vr~ x-^L • (9) 

I m fr — m r, ) 

In the following we will drop the subscript M when speaking of vrmi an d 
simply consider v R to have a small left-handed component when left-right 
equilibration processes are in equilibrium. Thus the RH sneutrino inherits the 
gauge and gaugino interactions of the LH sneutrino (see for example [18]): 

Xav 



£qauqino—Vfi — V^32 



m? — m? 



v c R v L W" + v%eW + + h.c. 



-V2g 17 — <-\i>%v L B + h.c] . (10) 



m? — m~ 



C g auge-9 R = ~i\ 7 ^ Mw ° ( W + U R * d g + W~? d V R 



1 1 

m~ — m~ 

-iX Z Zut d v R . (11) 

[mi L ~ mi R ) 

The list of the RH sneutrino interactions (excluding 4-points interactions) 
is included in fig.([T]). It is interesting to note that in the MSSM + Dirac 
(s)neutrino, while the RH neutrino remains completely sterile, the RH sneu- 
trino can mix with its LH counterpart and through it interact with the gauge 
sector. 

Before continuing we should mention what happens when one introduces 
a Majorana mass for the RH (s)neutrino. The superpotential is now 

W D XUeijHjN + M R NN (12) 

which results in the following added .F-terms in the Lagrangian: 

C F ,majorana = -2M 2 R D C R D R * + (M R Xeh + u c R * - M R XD L h°D c R * + h.c.) . (13) 

Adding the SUSY-breaking terms, 

£majorana= ~ {jn% R + 2M R ) V R V R - m| £ V* h i> L 

+ A (M R + a) (eh + D R * - 9 L h°9 c R * + h.c.) . (14) 

We want only to consider here the case of the RH sneutrino as the LSP; then 
automatically this excludes a very large, see-saw type of M R . A Majorana 
mass small enough for the RH sneutrino to be the LSP and that does not 
allow for thermalisation might be possible, as one can infer from the next 
section. Indeed it was shown recently that the right-handed sneutrinos can be 
non-thermal in the presence of a Majorana mass [20]. Such cases of pseudo- 
Dirac neutrinos are highly constrained from a phenomenological point of view 
(see for example [19]). We shall not consider this case any further. 
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|M| 2 = 8AVA 2 (p w>B ■ pi - m w ^m l/ ) \M\ 2 = X 2 M^ Z A 2 [ - (p u - R + pj) ■ (jp v ~ R + p~) 

+ vk~z ( Pv ~ R + ' P w < z (P"« + Pi) ' P w > z } 

Figure 1: RH sneutrino interactions with their corresponding amplitude, with A 2 
defined as A 2 = a 2 / (m 2 ^ — m |k) • 4-points interactions amplitudes are all simply 
proportional to A 4 



3 (Non-) thermalisation of the RH sneu- 
trino 

Let us first go back to the interactions of the sterile part of the RH sneutrino. 
Among these, only a four-point interaction or a LH sneutrino or a higgsino 
exchange contribute to RH sneutrino annihilation (see fig.([T|)). They will 
be out of equilibrium as long as their interaction rate T is smaller than the 
expansion rate H = T 2 /Mp. This necessitates for the 4-point interactions 
that 

r 4 ~ A 4 T < H => T > KT 33 GeV (15) 

where we've taken A = 10 -13 (considering the observed neutrino mass splitting 
squared to be indicative of the neutrino masses). This condition is obviously 
respected throughout the history of the Universe. The LH sneutrino exchange 
and higgsino exchange have rates Ty L ~ A 4 a ,T/rn~ L and T ^ ~ A 4 T 3 /m^.; 
imposing that these be smaller than H constrains the parameters such that 

To, < H => -V < 10 33 T 

m k 

r s < H => < lO^GeV" 1 . (16) 
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The first constraint is the strongest at low temperature while the second is 
strongest at high temperatures, so let us consider the first at the temperature 
no and the second at the highest reheating temperature that escapes the 
gravitino problem, T RH ~ 10 9 GeV [21-23]: 



a < KTm^GeV 



in 



ii 



> l(T 12 GeV . (17) 



The second condition is certainly met; as for the first one, it is a quite mild 
condition on the trilinear coupling a that can be easily met. Hence these 
interactions do not generally reach equilibrium. Moreover, before the elec- 
troweak phase transition, left-right mixing is out-of-equilibrium as well. These 
could go through 4-points interactions, which we already know to be out of 
equilibrium; they could also go through higgs exchange, which has a rate of 
^LR,h ~ \ 2 a 2 g 2 /T . Imposing the rate to be smaller than H all the way down 
to the electroweak transition, T ew forces 



o< \/AS;~ 107GeV ' <18) 

which again is easily respected. Hence the RH sneutrinos never equilibrate 
with the plasma before the electroweak phase transition. 

After the electroweak phase transition, we consider left-right mixing through 
mass insertiorJE The RH sneutrino now has an active part proportional to 
(Xav)/ {m 2 L — m 2 R )- There are now gauge and gaugino interactions that al- 
low for RH sneutrino annihilations; their rates are respectively given by 



a A \ A v A g 2 T 5 
m z ( m l L ~ m l R ) 



n 4\4 4 2 T 3 

% - a , XV9T N2 - (19) 



w V Ul u r 



Again, imposing these rates to be smaller than H is turned into a constraint 
on the unknown parameters 

< 10 6 GeV _1 . (20) 



2 2 

m~ — m~ 



This constraint implies that a mass degeneracy between left- and right-handed 
sneutrinos might eventually allow thermalisation of the RH sneutrino, de- 
pending on the size of the trilinear coupling. For a trilinear coupling and a 
LH sneutrino mass of order lOOGeV, only a RH sneutrino mass greater than 
~ lOOGeV - lOeV (but smaller than lOOGeV for it to be the LSP) can spoil 



2 Would left-right equilibration remain out of equilibrium after the electroweak phase transition, 
it would not allow for a quicker annihilation equilibrium and would only decrease the final relic 
density. Hence we will consider that after the electroweak phase transition gauge and gaugino 
channels become open as it is a safe assumption that will also simplify the calculation. 
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the non-equilibration; for the case of the mattergenesis mechanism of [16], 
a ~ lOOGeV and ~ lGeV, the LH sneutrino mass has to be between 
1 and 1,0001 GeV to allow thermalisation. The condition of eq. (|20p is thus 
fairly mild. 



< 3 x 10" 5 . (21) 



We can re-express eq, (|20|) as 

Xav 



2 2 

m~ — m~ 



This implies that the active part of the RH sneutrino (see eq.©) has to be 
small in order not to allow thermalisation. A LH part that's non-negligible 
when compared with the Yukawa coupling is still allowed however, and for 
this reason the gauge and gaugino couplings cannot be neglected altogether. 

So it appears that as long as the active part of the RH sneutrino is kept 
under control, the (LSP) RH sneutrino never thermalises. This is of little 
interest in a case where one is trying to simply have RH sneutrino DM with 
a relic density matching the observed one. It is however of interest in the 
opposite case of DM created through mattergenesis; supposing that a certain 
mechanism has created £lu R , we have shown up to now that this asymmetry 
will remain unaffected by the thermal plasma at later times. 



4 Dirac RH sneutrino relic density 

In the typical picture of a supersymmetric partner that decouples after becom- 
ing non-relativistic, the weak interactions of the RH sneutrino would generally 
produce a large relic density [24]. However the fact that RH sneutrino annihi- 
lations never equilibrate steers us away from the conventional calculation, and 
as we'll see the final relic density can be in fact very low. Here RH sneutrinos 
are created by the decays of the type 1^2 + vr, where particles 1 and 2 are 
in equilibrium. The relic density is given by solving the Boltzmann equation 
(see eg. [4]), which states here that the time evolution of the number of RH 
sneutrinos riy R follows 



+ mn„ R = Ci , (22) 



where the term linear in n„ R accounts for the expansion of the Universe, and 
Y^i Ci is the 'collision' operator. For the case at hand ^ Q is given by 

a = IU (pi -P2-P* R ) |m,i 2 (i ± h R ) (i ± h) h 

(23) 

where TL X = (2-7r) 3 2S x , f x are distribution functions (Fermi-Dirac or Bose- 
Einstein for particles 1 and 2), and there's a Ci to be calculated for every 
decay channel using the corresponding \M\ 2 as given on fig.([T]). Eq. (122ft can 
be rephrased in terms more convenient for a numerical solution using the yield 
variable Y = n/s, where s is the entropy density of the Universe. In terms of 
Y the Boltzmann equation becomes 

Y = I 7k¥ dT - < 24 > 

J- now 
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a 


m i>L 




m H 




100 


300 


100 


1000 


3 x 10" 


-5 


1000 


300 


100 


1000 


8 x 10" 


-3 


2000 


300 


1 


1000 


3 x 10- 


-3 


100 


300 


1 


1000 


3 x 10- 


-5 


100 


100 


1 


1000 


2 x 10- 


-4 


100 


50 


1 


1000 


2 x 10" 


-3 



Table 1: Various set of parameters and the relic density they generate. All masses 
are in GeV. The fourth line corresponds to the mattergenesis model of [16]; the 
scenario suggested in [9] only stated that a sneutrino mass of 0(1) GeV should be 
used (lines 3 to 7). The scenario presented in [17] requires a RH sneutrino mass 
of ~ lOOGeV. A large trilinear coupling or a close degeneracy of sneutrino masses 
enlarge the relic density. 



The yield variable is directly related to the density parameter through 

_ PvrI s _ m u R Y 

rehc ~ Q crit ~ 3.5 x 10-9/1 2 ( ' 

where £l cr it and h can be obtained using [25]. Going back to the collision 
operator, using the Maxwell-Boltzmann approximation [4] and integrating 
everything that can be integrated straightforwardly, we obtain 

l = ^f III E lE , R E 2 \M^ M N hS {El ~ E »* ~ E2) ■ (26) 



where Ei now stands for \j\p\ — Pv R \ 2 + m\. 

We have solved eq. ()24p for a number of sets of parameters that respect 
the constraints obtained in section [3] (tabled]). The relic densities obtained 
are completely negligible compared to the observed DM density. Using much 
larger trilinear couplings (but allowed by the non-thermalisation constraints) 
or degenerate sneutrino masses would lead to large relic densities [26] . It is not 
the case, then, that large relic densities are simply impossible; they can indeed 
be achieved by a certain level of parameter tuning. What we wish to suggest 
here, however, is that, as exemplified in table [lj for a large number of models 
that respect only mild constraints on the trilinear couplings and sneutrino 
masses, the relic density of RH sneutrino is indeed very small. When this 
observation is coupled to the fact that the RH sneutrino is not thermal, it 
leaves wide open the possibility of DM generation through mattergenesis. 

This is not the end of the story, however, because we have up to now only 
considered the case where all the other particles involved (apart from the RH 
sneutrino) are in equilibrium. They will eventually freeze-out, however, and 
since none of them is the LSP, their relic density will eventually be 'dumped' 
into the RH sneutrino one. In other words, the next-to-LSP (NLSP) (or 
MSSM-LSP) relic density should also be calculated and added to the RH 
sneutrino relic density [26] (see fig.©). This is the one thing small Yukawas 
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Log ( T ) ( GeV ) 

Figure 2: Evolution of the RH sneutrino relic density as a function of temperature 
(time running backwards). The parameters that have been used here are the ones 
in the fourth line of table [TJ The next-to-LSP will freeze-out at around typically 
at which point the RH sneutrino relic density has already reached its 
final value. Some time after the NLSP freeze-out the NLSP relic density will be 
'dumped' into a RH sneutrino one, thus adding a (hopefully small - see text) 'step' 
to this plot. 



give no protection against: dumping of large amounts of RH sneutrino by the 
decay of a MSSM-LSP that would happen to have a relic density comparable 
to the observed dark matter density. Fortunately a low MSSM-LSP is still an 
open possibility. Indeed it has been noticed recently that some LSP candidates 
usually considered to have 'correct' relic densities can have instead either 
much too high or - more importantly here - much too low relic densities (see 
for example [27] or [28] for an overview). It should be mentioned, however, 
that no matter what their DM candidate is, any mattergenesis scenario has to 
assess the question of MSSM-LSP dark matter, because if indeed the MSSM 
with R-parity is a reality, then the LSP is necessarily a source of dark matter. 

The left-right 'non-equilibration' before the electroweak phase transition 
was noticed in (non-supersymmetric) Dirac neutrinos in [29], where it was 
argued that it can be used to construct a leptogenesis mechanism in the neu- 
trino sector. The mechanism was used further in the MSSM in [16], where 
this time it rendered possible the creation of both baryonic and dark matter 
through the Affleck-Dine mechanism, a possibility also noticed more recently 
in [17]. The parameters that yield a correct DM density in [16] (a ~ lOOGeV 
and mjj R ~ lGeV, no constraint on m„ L ) are 'generic' enough to respect all 
the constraints of section [3] and produce a small relic density (see table [I]) . 
This ensures that the RH sneutrino density created around reheating time 
by the Affleck-Dine mechanism remains as such and can be straightforwardly 
taken as the DM density, as was actually assumed in [16]. It is interesting 
to note that while [16, 17] showed that the smallness of the Yukawa coupling 
provides a simple way of creating a DM density, what we've shown here is 
that it is also a built-in protection of this density. We also mention in table 
Q]the scenario of [9], in which a new, heavy Majorana fermion X decays to a 
RH sneutrino and a SM fermion. 
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5 Conclusion 



The Dirac RH sneutrino LSP is a natural candidate for DM within a matter- 
genesis scenario. Unless the trilinear coupling a or the degeneracy between 
left- and right-handed sneutrino become very large, the RH sneutrino never 
thermalises throughout the history of the Universe. This is a desirable charac- 
teristic in most mattergenesis scenarios as it ensures that a DM energy density 
created in the early Universe is not reprocessed at later times by fast anni- 
hilations. Moreover the relic density of RH sneutrino LSP is generally much 
smaller than the observed DM density, which implies that no large amount of 
additional DM is added to a density created earlier. Again this remains true 
as long as both trilinear coupling and sneutrino mass degeneracy are kept un- 
der reasonable control. In the mattergenesis scenarios suggested in previous 
works [9,16,17], the constraints obtained here are easily respected. As such 
the MSSM+z/r appears as the minimal extension to the MSSM that allows 
for DM to have mattergenesis as its source. 
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